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Multilayer "ferromagnet-layered antiferromagnet" 
(Fe/Cr) structures frustrated due to the roughness of layer 
interfaces are studied by numerical modeling methods. The 
"thickness-roughness" phase diagrams for the case of thin fer- 
romagnetic film on the surface of bulk antiferromagnet and 
for two ferromagnetic layers separated by an antiferromag- 
netic interlayer are obtained and the order parameter distri- 
butions for all phases are found. The phase transitions nature 
in such systems is considered. The range of applicability for 
the "magnetic proximity model" proposed by Slonczewski is 
evaluated. 
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I. INTRODUCTION 

Multilayer magnetic structures have attracted a 
widespread attention after discovery of the phenomenon 
of giant magnetoresistance.El Within the span of eleven 
years a great many of papers devoted to the multilayer 
structures has been pubhshed. Considerable recent at- 
tention has been focused on multilayer "ferromagnet- 
layered antiferromagnet" structures. An example of such 
structures, according to recent experiments on ncuU-on 
diffraction, is given by multilayer Fe/Cr structures ,BH in 
which antifcrromagnetic chromium layers of the thickness 
d < A.bnm consist of ferromagnetic atomic planes with 
antiparallel orientation of spins in adjacent planes. The 
spins of chromium atoms lie in these planes, the plains 
in their turn being parallel, on the average, to the layer 
interfaces. An analogous structure has been obsecKed in 
manganese layers in multilayer Fe/Mn structures Jj'El 

In multilayer "ferromagnet-layered antiferromagnet" 
structures, the exchange interaction between ferromag- 
netic layers is caused mainly by the interaction via anti- 
ferromagnetic order parameter, whereas the RKKY inter- 
action far away from the Neel temperature provides only 
a small additive to the first oneJs To describe the interac- 
tion via antifcrromagnetic order parameter, Slonczewski 
proposed a phenomcnological model of "magnetic prox- 
imity" in whose context the layers of ferromagnet are 
assumed to be magnetized practically uniformly and con- 
siderable distortions of the order parameter are supgosed 
to take place only in the layers of antiferromagnet. Q 



The presence of atomic steps at layer interfaces that 
change the local thickness of the antiferromagnet by one 
monoatomic layer gives rise to frustrations in the system 
(Fig. 0a). In this case the uniform distribution of the 
order parameters in the layers does no longer correspond 
to the energy minimum. 

If the characteristic distance between atomic steps at 
the layer interface (we shall call it the step width R) 
exceeds some critical value, it becomes enejcsetically fa- 
vorable to break the layers up into domainsErB (Fig. |l|b). 
The domain boundaries in layer planes coincide with the 
atomic step edges. It should be noted that the magniiiide 
of R strongly depends on the technology conditionsEj. 

Recent investigations on the state of ferromagnetic iron 
film on the rough (001) chromium surface have shown the 
presence of several magnetic phases depending upcp, the 
film thickness and roughness rate (the value of R) .llil 

All the aforesaid bears witness to considering the 
problem of "thickness-roughness" phase diagram of 
"ferromagnet-layered antiferromagnet" magnetic struc- 
tures as a high priority task. 

The paper has the following structure. The consid- 
eration begins with a simple model that allows one to 
describe the behavior of the system in question qual- 
itatively. In the third section one can find model re- 
sults on the phase diagram of the two-layer system: 
thin ferromagnetic film on antifcrromagnetic substrate 
(or vice versa). The phase diagram for "ferromagnet- 
antiferromagnet-ferromagnet" three-layer structure is 
found and discussed in the fourth section, the treatment 
being possible to be generalized to multilayer. In conclu- 
sion the main results are adduced. 



II. DESCRIPTION OF THE MODEL 

When describing the multilayer structure, we restrict 
ourselves to the mean-field approximation. Let us in- 
troduce the order parameter for each magnetic layer: 
the magnetization vector for the layers of ferromagnet 
and antiferromagnetism vector, equal to the difference in 
magnetization vectors of the sublattices, for the layers of 
antiferromagnet . 

It is known that in sufficiently thin magnetic layers 
with the thickness of some nanometers, the atomic spins 
lie in the layer planes. Therefore, at T < Tq, Tm (here 
Tc is the Curie temperature for single ferromagnetic lay- 



1 



ers, and T/v is the Neel temperature of antiferromagnetic 
layers), it is possible to describe the local value of the 
order parameter in the layer planes by the angle of the 
order parameter vector with a certain given axis in the 
layer plane. The order parameter modulus is assumed to 
remain unchanged in each layer. 

In the framework of the assumptions made above, the 
exchange energy Wi caused by a nonuniformity within 
the i— th layer can be introduced as 

where the integral is taken over the i— th layer volume, 
Ji is the exchange stiffness, Si is the mean spin of the 
atom, bi is the inter-atomic distance. 

By varying Eq. (|^) with respect to 6i we get the equa- 
tion for the order parameter distribution in the layer: 

Ae, = 0. (2) 

To obtain the true boundary conditions one needs a more 
exhaustive procedure. Owing to frustrations arising at 
the layers boundaries the difference 6i — 9i^i can be- 
come large enough at the interlayer boundary, whereas 
inside the layers the frustrations are absent, the value 
of 9i varies smoothly, and the difference of 9i values in 
the nearest cells is small. That is why, when calculating 
exchange energy within the X-Y model, we may expand 
cosine of the difference of the angles 9i in neighboring 
cells in power series if the cells belong to the same layer, 
but we'll not do that if the cells belong to different layers. 
By differentiating the energy with respect to the quantity 
9i in the cell lying in the layer boundary, we come to the 
equation which takes the following form as one passes to 
the continuous representation: 

where A is the two-dimensional Laplacian in the layer 
plane, ^ is the derivative in the direction of the outer 
normal to the interface plane, the exchange constant 
Jf,af describes the interaction of spins corresponding to 
different layers, and all distances are normalized to the 
value of bi = b which is assumed to be the same in all 
layers. The signs in the right-hand part of Eq. (H) arc 
opposite for different sides of the atomic step at the in- 
terface. For the free surface of the layer the right-hand 
side of Eq. (||) equals zero. 

The energy of exchange interaction between adjacent 
layers takes the from 

T^M+i = / cos {9, - dS, (4) 

where the numbers i and i + I denote corresponding lay- 
ers, the integration is performed over the interface of the 
layers, and the sign in the right member coincides with 
that one in Eq. (||). 



By varying the energy of interlayer interaction repre- 
sented in terms of the continuum model distribution one 
gets an equation similar to Eq. (||) but without the first 
summand in the left-hand side. This renders it impossi- 
ble to pass from Eq. (||) to Eq. (||) if adjacent layers are 
identical in composition. 

The atomic steps divide the whole interlayer surface 
into regions of two types: in the regions of the first type 
the surface energy takes its minimum for 9i = 9i+i and in 
the regions of the second type the energy is a minimum 
for 9, ^ TT - 6'j+i, 

To find the order parameters distributions in the mul- 
tilayer structure it is necessary to solve the system of 
differential equations (||) with boundary conditions (||). 

Now let us take up the applicability of this simple 
model to real multilayer structures. The continuum ap- 
proximation is valid only in the case that characteristic 
measures of the problem are much higher than the in- 
teratomic distance. The layer thicknesses in multilayer 
structures as well as the distances between atomic steps 
are of about several nanometers, so they may be consid- 
ered as far exceeding the lattice parameter. Hence the 
continuum approximation is applicable for estimations 
and qualitative treatment. 

In the model proposed, the exchange interaction is as- 
sumed to be isotropic, i.e. identical in the layer plane and 
in perpendicular direction. The model with anisotropic 
exchange interaction is reduced to the given one by 
renormalization of the distance scale in any of two non- 
equivalent directions. 

An interdiffusion of atoms from neighboring layers 
leads only to renormalization of J/,a/, if the region of 
mutual stirring extends up to several monoatomic lay- 
ers, i.e. is of atomic size scale. Thi|S-.constant value is 
obtained by microscopic calculations £3 

Eqs. (||), (H) are written in the exchange approxima- 
tion, but can be easily generalized to the case of small 
anisotropy in the layer plane. 

III. THE METHOD OF CALCULATIONS 

Let the atomic step edges be rectilinear and parallel 
to each other. The a;-axis of the coordinate system lies 
in the layer plane and is perpendicular to the step edges, 
and the z-axis is perpendicular to the layer plane (two- 
dimensional case). 

The assumed set of equations includes Laplace equa- 
tions (^) for the each plane layer — oo<a;<oo, < z < 
Ui, where i — 1 ... n is the layer number and is the layer 
width, the boundary conditions (^) at interlayer surfaces 
being nonlinear and discontinuous. For numerical calcu- 
lations we reduce this set to the system of unidimensional 
integral equations. 

We require the function 9i{x, z) to be continuous inside 
the region 0<a;<L, Q < z < ai (here L is half of the 
modeling region length along x axis), and additionally 
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^L_>OL ~^ ^' '^'^ make the function 9i{x,z) periodic 
with the period 2L we also put 6i{—x, z) — 9i(x, z), —L < 
X < L. Let us put a uniform net {xj} with cells number 
N and step = 2L/N on the interval —L<x<L and 
introduce the function 9i{x, z), as well as the right-hand 
part of Eq. (S) as the Fourier series: 



N 



,{x,z) = 4>o,j(2;) + 2^*fe,j(2)e*i 



fc=i 



a±(x)sin(5^±(a;))=vl/± + 



N 
k=l 



(5) 
(6) 



Here (x) = ± 



Jf.ajSi±l 
JiSi 



is a stepwise function changing 



its value at the step boundaries and 

60, = 69+ (x) = -6e-^,{x) = 9,+iU=o 



Substituting Eqs. (|5[]6|) into the starting equations one 
comes to normal differential equations for any harmonic 
k = 1, 2 ... TV inside each layer i: 



dz^ 



k^k 



0, 



(7) 



(the index i is omitted here and further, where it does not 
make understanding difficult) with the following bound- 
ary conditions 



dz 



± 

k ' 



(8) 



z—a.O 
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A, = -^l-cos(|fcA^). 
The solution of Eq. (|^) is sought in the form 

with the constants Cf , derived from boundary con- 
ditions (ph. As a result one has 



Mz) = K+{z)^+ + K^iz)^^, 



(9) 



where 



K+{z) = dKke-^--^'^^- [(1 + Ak) + Ak)e 



2zAk 



K^{z)^dKke 



(l + Afe) + (l-Afe)e-2(---)^^ 



dKk 



-A^'[{l + Akf-{1-A,fe 



The equation for zeroth-order Fourier coefhcients has the 
form 



= 



dz^ 



with boundary conditions 
dz 



With account of = for free boundaries of the top 
and bottom layers, and also of the constant ratio 

for the neighboring layers, one has = 0. If one sepa- 
rates out the average angle 

39 = ^ f 69{x)dx 
2^ J -L 

and the variation 

69 {x) = 69{x) - 69, 
it is possible to find 

= / a{x)sm{69{x) + 69)dx 



cos 69 I a{x) sin 69 {x)dx 
sin 69 I a{x) cos 69 {x)dx. 



Hence the average angle is 



69 — jm — arctan 



/_L '^{^) sin 69{x)dx ^ 
^ (j{x) cos 69{x)dx ^ 



(10) 



With Eqs (^|6|) and (|9[]10|) one can obtain the desired 
integral equation for the function 69i: 



N 



1 



k=l 



e-'i'^^d^ 



2L j_L 

K^+iJ^W+i + K++i,ki^H+i sin 69,+i 

~Kr',^{ai)al sin69i^i - K^j^{ai)al sin69i 

= i{K, 69„69,±i). 

To solve this equation numerically we use a simple itera- 
tion procedure 



69, 



,n+l 



il-F,{x)) 69': + F,{x) I{Fk^,K, 



with < F(a;), Ffc < 1 which are filters providing stabil- 
ity of the procedure and increasing its convergence rate. 
The iteration procedure is over when the residual solution 



max 



\x)-I{69-,69-^,) 



becomes less than a given value (normally, e < 10~^). 

Then the solution over the whole region can be recon- 
structed from the formula: 



i{x,z) 



z—a,0 



L 

k=l —^-L 

Kik{z)<yi svn59i_i + K+f^{z)(7+ sin69. 
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As an initial approximation, one can take the function 
of the form 



where r]{x — Xj) is a unit stepwise function with a jump 
at the position of crystal lattice defect location xj and 
the sign ± means that the given function can enter the 
sum with both plus and minus signs. As a result, one 
has ^ possible initial approximations, where Nj is 
the number of defects in the region under calculations. It 
should be noted that to each of the initial conditions there 
can correspond one of the solutions of initial nonlinear 
equation (a local minimum of potential energy). To find 
the global minimum, it should be compared the energies 
corresponding to all of the solutions obtained. 

The solution depends on the values of a^, characteris- 
tic distance R between the edges of the steps, and the 
parameter 



Jf,afSaf 



(11) 



characterizing the ratio of the exchange interaction ener- 
gies of the nearest spins belonging to different layers and 
to ferromagnetic layer, respectively, and also the param- 
eter aaf defined by Eq.(|ll|) with the subscript substitu- 
tion: / <-> af. 



IV. THIN FERROMAGNET FILM ON 
ANTIFERROMAGNET 

Let us consider a thin ferromagnet film of the thickness 
a (in normalized unitless scale) on the surface of layered 
antiferromagnet. The problem of antiferromagnet film 
on ferromagnetic substrate is easily reduced to the given 
one. 

Three different states of the film-substrate system are 
possible, depending on the relationship between indi- 
cated parameters. If the atomic steps on the film- 
substrate interface are well away one fromj-ajiother, the 
frustrations give rise to domain formation.Ercl The mag- 
netization orientation in each domain is conditioned by 
the surface energy minimum (Fig. |b). Hence the mag- 
netizations in the neighboring domains are antiparallel. 
The structure of Neel walls demarcating the domains de- 
pends on the value of dimensionless parameter 7 equal 
to the ratio of exchange energies in the film and in the 
substrate: 



7 ^ 



JfS' 



f 



af ' 



(12) 



We have already considered the case where the orjier pa- 
rameter distortions in the substrate are neglig ible.ElSuch 
an approximation conforms to 7 <C 1. A distinguishing 
feature of the domain wall arising is an increase in its 



width with a distance from the edge of the atomic step 
giving rise to the wall. The parameter afU therewith 
plays an important part in the wall width behavior. The 
width of the domain wall S(z) is interpreted as the dis- 
tance between the points with coordinates (xi,z) and 
(x2,z), which correspond to 61 = 7r/4 and O2 — 37r/4, 
respectively. 

If a/ a <^ 1, the wall width S dependence on z is of 
no significance and one has a one-dimensional problem. 
The wall width magnitude Sf can be estimated from the 
following simple consideration. 

The \V9\ value inside the wall is of the order of Sj^. 
Estimating the energy Wi (Eq. ^ per the domain wall 
unit length, one obtains 



Wl 



(13) 



At the same time, the spins on the interface are frustrated 
in the region |x| ^ Sf, resulting in the energy increase in 
excess of the minimum by the value 



W2 



Jf,afSfSaf5f 



per the domain wall unit length. 

Minimizing the sum w = Wi + W2 one finds 



\/a/a/: 



and the domain wall energy per unit length is 



— ^JaJfJf^afSfSaf. 



(14) 



(15) 



(16) 



The precise numerical calculations of 5/ and w values 
for a wide range of parameters af and a confirm the 
estimations given above (the same relates to all cases 
considered below). 

If a/ a 3> 1, the domain wall width essentially increases 
as compared to the value of 5^ at the interface. The 
5f{z) characteristic dependence is shown in Fig. ^. The 
increase is practically linear [85 /dz « 1), with slowing 
down nearly to zero at a free surface of the film. 

One can find the values of 5^ and w using the approach 
analogous to that proposed for the case a/a ^ 1. In the 
region 8^ <^ p <^ a the value of jV^j is proportional 
to p^^ , where p is the distance from the step in the izlz 
plane, just similar to the Kosterlitz-Thouless vortexo. 
The value wi is equal to 



Wl 



2 b 



In 



(a/5i 



(17) 



and 'W2 is given by Eq. |l^ with substitution of 5^ instead 
oi 5f. After minimization one has 



af 



af 

5f{a/2)^a, 



(18) 
(19) 
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and w K wi ^ W2- 

The value of is of the order of interatomic distance, 
and a mean wall width comprises tens of angstroms, i.e. 
the domain walls resulting from frustrations are much 
more narrow than "normal" domain walls in ferromag- 
ncts whose width is due to a competition between the 
exchange and anisotropy energies. 

In the case of iron film on chromium substrate consid- 
ered in this paper, the value of 7 is high, 7 3> 1. For such 
conditions the domain wall structure is more complicated 
because the order parameter distortions can extend to the 
substrate as well. 

If aafja <C 1, the order parameter distortions in the 
substrate are small and the domain wall characteristics 
do not differ from those in the case 7^ 1, a/a 1. 

For aafja ^ 1, as one can see from the results of 
modeling (Fig. ^), two characteristic lengths come into 
account: the first one 6q^ is the domain wall width in 
the substrate near the interface, and the second one Sf is 
the domain wall width in the ferromagnetic layer. Since 
(5/ S> a, one can neglect the domain wall widening in the 
ferromagnet. 

Let us consider the order parameter behavior in the 
antiferromagnet and the situation at the film-surface in- 
terface. Let the atomic step on the interface coincide with 
the y-axis of Cartesian coordinate system. For x <C —5/ 
the conditions 6a f = Of = are fulfilled, and for x ^ 5/ 
the conditions Oaf = 0, Of = tt are met. It follows from 
the symmetry of the problem that Oaf =0, Of = tt/2 for 
X = 0. The width of the region at the film-surface inter- 
face where the value of ^/ — Oaf differs from its optimum 
(0 for a; < and tt for x > 0) equals Sq-^ . In the region of 
I a: I ^ 6f and \z\ < 6f there arise distortions of the order 
parameter in the substrate (Fig. ^). 

Analogously to the previous consideration one can esti- 
mate the energy wi inside the film by Eq. ( p^ and inside 
the substrate as 



af ^ ^ '^"■f^af , Sf 
~ — ; ^In^^r, 



(20) 



The energy W2 is given by Eq. (|lj) with Sq^ instead of 
Sf. After minimization one has 



.a/ _ l + agf 

On ~ 1 



aaf 



Sf « 7a, 



(21) 
(22) 



and w « w'^'^ , so the main contribution being caused by 
the order parameter distortions in the substrate. 

The foregoing estimations relate to the case of well 
away distances between the steps. If the distance be- 
tween the steps decreases and becomes less than the crit- 
ical value Rc — Sf, the domain walls begin to overlap and 
the film switches to a single domain state. The transition 
from a polydomain state to a single-domain state is con- 
tinuous and, in the strict sense, is not a phase transition. 



If 7 <C 1, a/a » 1 and (5g <C i? <C <5/, a static spin 
vortex arises in the film near the substrate enclosing the 
area z ^ R. 

But if 7 3> 1, aafaj 3> 1, and the value of R is in the 

range Sq'^ <C R ^ Sf, a similar spin vortex arises near 
the interface in the substrate. Near the interface spin 
orientation in the vortex corresponds to the minimum 
of the interface energy W2- Far from the interface spin 
orientation is homogeneous. The results of modeling are 
introduced in Fig. 

For all others cases the distortions of order parameters 
are small in both the film and the substrate and their 
values can be considered as constants. 

As it was already mentioned, the steps divide the whole 
interface into the areas of two types. If the mean mag- 
netization vector makes an angle ip with the antiferro- 
magnetic order parameter in the substrate volume, then 
the value oi Of (or Oaf) changes from zero to ip in the 
vortex occupying the first type area, whereas in the vor- 
tices occupying the second type areas the value of Of (Oaf) 
changes from ip to n. j-. 

By an analogy to the "magnetic proximity" modefl the 
energy of the system can be written down as 



(23) 



where the ratio of phenomenological constants Ci and 
C2 is proportional to the ratio of the areas of both types 
((Ti and (72 correspondingly) . If the distributions of these 
areas in sizes for each type are equal, one has 



Rb 



(24) 



In the case of cti = (T2, without regard for the energy 
of anisotropy induced by the steps, in both the vortex 
phase and the region of small distortions the equilibrium 
magnetization of the film must be perpendicular to the 
antiferromagnetic order parameter. The phase diagram 
of the two-layer system is exhibited in Fig. ^. 

This specific pattern correlates well with the results 
of paperjlil where the "thickness-vicinal angle /?" phase 
diagram for iron film on Cr (001) surface was investi- 
gated. For /3 angles close to zero a polydomain phase 
was observed for the film thickness a < ac = 3.5nm. For 
the film of the critical thickness Uc the characteristic dis- 
tance between the edges of random steps fits the value 
-fac- In thicker films a single-domain phase was observed, 
with magnetization perpendicular to the step edges. The 
theory proposed above predicts the antiferromagnetism 
vector to be parallel to the step edges, its experimental 
verification being of immediate interest. 

If /3 7^ 0, regularly distributed parallel steps are added 
to random atomic steps. When the concentration of 
regular steps becomes dominating (3 > 1°, the critical 
thickness Oc falls down. It follows from our theory that 
ttc cc R (X tg~^(3 (X 

At high values of P the orientational phase transition to 
the phase with magnetization parallel to the steps takes 
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place.0 This transition is caused by the anisotropy in- 
duced by the steps thrp«gh relativistic effects, for exam- 
ple, dipole interactionJla 



V. MULTILAYER "FERROMAGNET - 
ANTIFERROMAGNET" STRUCTURE 

Let us consider for simplicity a three-layer system con- 
sisting of two ferromagnetic layers and an antiferromag- 
netic interlayer. In view of a number of various param- 
eters, we restrict ourselves to the case of 7 ^ 1 and the 
layers of equal thicknesses. Such a system can occur in 
three different states. 



^-phase 



this region |V6'| cx a ^. The contribution of the region 



to W, ■' IS 



(2) _ J<^f^lf^'f 

'1 ^ ba ■ 



(26) 



Taking into account the contribution to wi from ferro- 
magnetic layers (Eq. (|l^)) and the interface energy W2 
one comes to the result given by Eq. (|l]) for 5^-^, 



and 



5'f ~ aV7, 



JafSlj ( a 



(27) 



(28) 



At high values of R all layers are divided into domains 
with parallel and antiparallel orientation of magnetiza- 
tions in ferromagnetic layers. The domain walls pene- 
trate into all three layers, their coordinates in the layer 
plane coincide with the atomic steps edges on each of 
two interfaces. The magnetizations of ferromagnetic lay- 
ers rotate to different sides in the domain wall The anti- 
ferromagnetic order parameter rotates together with the 
magnetization of that ferromagnetic layer whose interface 
is free of a step in the given place. 

The structure and energy of the domain wall depend 
on the parameter aafa. When aafa <^ 1, one can neglect 
the dependence 9f(^af){z), i.e. the domain wall widening 
(Fig.|a). 

The consideration analogous to the case of the film 
gives the expression (15) for the domain wall width in 
the ferromagnet and Saf ~ \J o-lotaf = 5fl\/l ^ ^/ for 
antiferromagnet one. The main contribution to the do- 
main wall energy wi comes from the ferromagnetic layers 
and w value is given by Eq. (|l6|). 

But if aafa ^ 1, then, as in the case of a two-layer 

system, two length scales come into existence: Sq'^ = 
(1 -I- aaf)/aaf and S'j: which differs from that given by 
Eq. (^H) because the distortions of the antiferromagnetic 
order parameter are now limited by the layer thickness. 
The distribution of the order parameters in this wall is 
shown in Fig. |^. 

The domain walls in ferromagnetic layers have the 
width much higher than a, thus it may be considered 
as a constant across the layer. In antiferromagnetic in- 
terlayer one can see two regions (Fig. ^) . In the first one 
a the situation is analogous to the two-layer sys- 

from the region 



X 



< 



tern: |V0| cx p ^. The contribution w 
is 

(1) ^ TT JafSlf, a 



w 



In- 



(25) 



But in the second region a <C |a;| <C (5j the Hues of con- 
stant values are almost parallel to the interfaces. In 



B-phase 

Since the magnetic rigidity of ferromagnetic layers is 
higher than that of antiferromagnetic (7 3> 1), the tran- 
sition to a state with practically uniform ferromagnetic 
layers, because of the domain walls overlapping, takes 
place with the distance R decreasing and reaching the 
value Rc = 5f{5'j). An additional energy is associ- 
ated with either distortions in antiferromagnetic inter- 
layer or peculiarities of the structure at the interfaces. 
Near the Neel temperature of an interlayer Tm {Tn is 
less than the Curie temperature of ferromagnet) one has 
7 (x Tjs! /{Tjs! — T), so the transition A ^ B may be per- 
formed by heating the system from an initial temperature 
To < T/v. Tie "magnetic proximity" model proposed by 
SlonczewsklD is appropriate just in that R range where 
the -B-phase exists. 

In the range of distances max ^a, 5^^^ ^ i? <C Rc, the 

energy dependence on the angle ■0 between the magneti- 
zations of ferromagnetic layers is described by Eq. ([2^). 
Really, for aafa ^ 1, in the regions of the first type 
(where parallel orientation of ferromagnetic layers mag- 
netizations is energetically favorable) the antiferromag- 
netic order parameter changes linearly with z from one 
interface to another by the t/i value (|V0| = ip/a). In the 
regions of the second type (where antiparallel orientation 
of ferromagnetic layers magnetizations is energetically fa- 
vorable) it changes hy it — ip, |V6'| = (tt — '4')/a. The 
frustration at interfaces is absent. The boundary energy 
between regions of the first and the second types is of no 
significance in the case i? 3> a. 

In the case aafa <^ 1 the antiferromagnetic order pa- 
rameter is almost constant, but one has frustrations at 
the interfaces. The jump of 6 value equals ip/2 at each 
interface in the regions of the first type and (tt — ip)/2 in 
the regions of the second type. The excess energy W2 per 
two-dimensional unit cell in the regions of the first and 
the second type equals: 
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W = , „ ' ^ sm — , 



62 



(29) 
(30) 



Eqs. ( 129^301 ) are not identical to Eq. (|23|), the phenomeno- 
logical constants Ci , C2 can be estimated by compar- 
ing energy difference of collinear and 90° orientations in 
the model proposed here and in the magnetic-proximity 
modeB: 



C1.2 



Cl 2 



tta/a > 1, 



2a 

4(\/2 — 1) 7 c c 2 „ 1 

^-p-^ Jf,af!SfSaf -p-, aa/a < 1, 



(31) 



where cri,2 are the areas of the first and second types on 
the layer surface. 

For (Ti = (72 the energy minimum is attained at ■0 = , 
i.e. the magnetizations of ferromagnetic layers are mutu- 
ally perpendicular in the absence of an external magnetic 
field. 

Ferromagnetic layers are, to-be sure, not ideally homo- 
geneous. By analogy to Ref.Ej the change of the quan- 
tities 9i, ip across the length of R can be shown to be of 

the order of [i?/(5/(5})j < 1. 

Since the radius of the vicinity influencing the value of 
9i at a given point equals (5/ ((5^), the deviation of 0i and 
-0 from their mean values is defined by the fluctuation of 
the quantity that characterizes the relation between the 
square areas of the regions of two types. If the regions of 
the first type dominate in given vicinity, one has ip < 7r/2, 
but if those of the second type, one has ^ > Tr/2. 

The characteristic area bounded on by the steps on 
the layer surface is of the order of R^, so the number 



N of the regions in the area Sf{6'f) 



is of the order of 



Sf{S'j:)/R . When it is assumed that the regions of both 

types are randomly distributed the predominance of the 
regions of one- type measures ViV « 5f{5'^)/R. Therefore 
the characteristic fiuctuation of the quantities 9i and 
is of the order of R/5f{5'f) < 1. 

In the ^-phase the energy of the system does not de- 
pend on the direction of the order parameter rotation in 
the domain wall. A different picture is observed in the 
-B-phase: when the walls overlap one another, the de- 
generation in the rotation direction is eliminated and a 
number of metastable states arise. They differ from each 
other by the direction and the value of the angle of rota- 
tion of antiferromagnetic order parameter in individual 
areas bounded by the atomic steps. 

And what is going on with further R decreasing? If 
Oiafa 1 and R lies in the range a <C i? <C 5a f, the 
system arrives at the region of small distortions where 
the order parameters are practically uniform, the mag- 
netizations of ferromagnetic layers being mutually per- 
pendicular and the constants Ci^2 lowering by a factor of 
{R/Saf)'^ comparing to their estimations (|3 



C-phase 

Now let us consider the range of small R, R ^ a. For 
such a case all distortions are localized near the inter- 
faces, the interaction between ferromagnetic layers be- 
comes weak and the main part is played by the energy of 
interaction between the neighboring layers demonstrated 
above with an example of a two-layer system. As the 
result, at cri = cr2 the antiferromagnetic order parame- 
ter is oriented perpendicularly to the magnetizations of 
ferromagnetic layers which thus appear to be collinear 
(parallel or antiparallel). Such a state is termed here as 
C-phase. 

For aa/a 3> 1 and the distance Sq^ <C i? ^ a, the 
static vortices are formed in an antiferromagnetic inter- 
layer near the interfaces. For lower values of R the system 
changes over to the region of small distortions. 

If aafCL <^ 1, the transition from the -B-phase to C- 
phase takes place already in the region of small distor- 
tions. Along with the -B-phase, the C-phase is character- 
ized by a presence of a wealth of metastable states. The 
numerical calculations has shown the i? — > C phase tran- 
sitions to be of the first order. Both phases coexist in the 
entire region of R values and their specific energies be- 
come equal at a certain value R* ^ a (Fig. The phase 
diagram of the three-layer system is exposed in Fig. ||. 
InliZI the distribution of spins in a two-layer ferromagnet 
— layered antiferromagnet system is modeled numeri- 
cally on the bases of Ising model. However, the Ising 
model corresponds to very strong anisotropy of easy-axis 
type and is invalid for the description of multilayer struc- 
tures of the Fe/Cr type, for which the anisotropy energy 
in the plane of the layers is much lower than the exchange 
interaction energy. The domain walls formed within the 
framework of the Ising model have atomic widths, so that 
the unique properties of domain walls in the multilayers 
ferromagnet — layered antiferromagnet where not ob- 
served in the cited paper. The orientation of spms in the 
three-layer system has been calculated in Ref.ll3 within 
the framework of a discrete model for the case i? ~ a, the 
magnetizations of ferromagnetic layers being assumed to 
be mutually orthogonal. The behavior of the system in 
the whole range of i?-values has not been analyzed. 

The value of R* is independent of temperature, there- 
fore the B ^ C phase transition is not observable un- 
der changes of temperature of the system. The tran- 
sition from the state with strong biquadratic coupling 
to the low-tempjesature state with a weak coupling be- 
tween the layeralH is unrelated to the transition consid- 
ered above. The transition considered in Refo is caused 
by the fact that as approaching T/v, the interaction of 
ferromagnetic layers via the antiferromagnetic order pa- 
rameter decreases and becomes equal to the interaction 
via spin polarizatioiuinduced in the antiferromagnet (the 
RKKY interaction) .H For higher temperature the latter 
dominates. 
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VI. CONCLUSIONS 

Let us state the main results and conclusions of the 
work exposed above. 

1. It is proposed a simple model allowing one to 
find the distribution of spins in "ferromagnet- 
antiferromagnet" frustrated layered structures. 

2. The "thickness-roughness" phase diagrams for a 
thin ferromagnetic film on antiferromagnetic; sub- 
strate and for a "ferromagnet-antiferromagnet- 
ferromagnet" three-layer system are obtained. 

3. The transition from a polydomain state to a single- 
domain state is continuous and, in the strict sense, 

is not a phase transition. 

4. The transition from the phase with mutually per- 
pendicular orientation of magnetizations of fer- 
romagnetic layers (_B-phase) to the phase with 
collinear magnetizations of the layers (C-phase) 
constitutes a true first-order phase transition. 

5. The "magnetic proximity" phenomenologi- 
cal model proposed by Slonczewski is shown to be 
adequate for the S-phase only. The parameters of 
the model are calculated for the entire range of its 
applicability. 
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FIG. 4. Static spin vortices in the film-substrate system for 
7^1. Levels of order parameters 9i in radians are shown 
by various hatching (see insert). The distribution was found 
for 7 = 8, a = 0.01, a = 8. The value 2 = corresponds to 
the interface between the film and the substrate. Steps are 
situated at the points x = ±10. 

FIG. 5. Phase diagram of the two-layer system. For 
demonstration the lines corresponding to equations R = 6/ 
and R = Sn'^ are drawn. 



FIG. 6. Domain wall in the three-layer system: Qa/a <C 1 
(a); Oaf a S> 1 (b). Levels of order parameters 9i in radians 
are shown by various hatching (see insert). The values z = 
and 2; = 16 correspond to the interlayer interfaces. Step is 
situated at the point x = 0, z = 0. The distributions were 
found for 7 = 10, a = 0.01, a = 16 (a) and 7 = 10, a = 1, 
a = 16 (b). In the Fig. (c) is represented the central part of 
the distribution shown in the Fig. (b). 
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FIG. 7. Dependence of B- and C-phases specific energy on 
the distance R between the steps (o = 64, a/ = 8, aaf = 1)- 

FIG. 8. Phase diagram of the three-layer system. For 
demonstration the fines corresponding to equations R = 5f, 
R = a, and R = Sq'^ are drawn. 
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